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Abstract. Josephson junctions and junction arrays are well studied devices in superconduc- 
tivity. With external magnetic fields one can modulate the phase in a long junction and create 
traveling, solitonic waves of magnetic flux, called fluxons. Today, it is also possible to device 
two different types of junctions: depending on the sign of the critical current density jc ^ 0, 
they are called 0- or tt -junction. In turn, a O-tt junction is formed by joining two of such junc- 
tions. As a result, one obtains a pinned Josephson vortex of fractional magnetic flux, at the 
O-TT boundary. Here, we analyze this arrangement of superconducting junctions in the context 
of an atomic bosonic quantum gas, where two-state atoms in a double well trap are coupled in 
an analogous fashion. There, an all-optical O-tt Josephson junction is created by the phase of a 
complex valued Rabi-frequency and we a derive a discrete four-mode model for this situation, 
which qualitatively resembles a semifluxon. 
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1. Introduction 

During the past two decades, the field of cold atomic gases has come a long way starting 
from almost lossless trapping and cooling techniques ^ to reaching quantum degeneracy 
of Bosons and Fermions |l2l. Many phenomena that are the hallmarks of condensed matter 
physics, whether in superfluid or superconducting materials jSl, are revisited within this novel 
context. Due to the remarkable ease with which it is possible to isolate the key mechanisms 
from rogue processes, one can clearly identify phase transitions, for example, Bose-Einstein 
condensation, the Mott phase transition or the BEC-BCS crossover. Today however, 
degenerate gases are still at a disadvantage if we consider robustness, portability or the ability 
for a mass production compared to solid-state devices, which is the great achievement of the 
semiconductor industry. Strong attempts to miniaturize cold gas experiments flU [51 and to 
make them portable [IH |71 [8l HI are currently under way in many laboratories. 

Due to the great importance and practical relevance of the Josephson effect in 
superconducting systems [[T0l[m[T2| . it has also received immediate attention after the first 
realization of Bose-Einstein condensates [[Il[Il[l5l[ia[l7l[Il[Il[20l|2B[2l[23IMl2l^ 
In particular, the combination of optical lattices with ultracold gases Il26ll27l has boosted the 
possibilities to investigate junction arrays experimentally. Remarkably, even the absence of 
phase-coherence between neighboring sites can lead to interference as demonstrated in ll28l . 
The possibility to study atomic Josephson vortices in the mean field description was raised 
first in connection with the sine-Gordon equation ||29l[30ll . 

In the present article we will report on such a transfer of concepts from a superconducting 
device [[3T1l . i.e., in various realizations of Josephson junction arrays and their unusual 
state properties of traveling (fluxons) and pinned (semifluxons) magnetic flux quanta to an 
analogous set up for neutral bosonic atoms in a trap. In particular, we will investigate an all- 
optical O-TT Josephson junction that can be created with a jumping phase of an optical laser. 

This article organized as follows: in Sec. [2l we will give a brief review of the current 
status of the superconductor physics of Josephson junctions. In particular, we will refer to the 
most relevant publications in this thriving field of fluxon and semifluxon physics; in Sec. [3l 
we will discuss a similar setup, which allows to find a pinned semifluxon in an atomic O-vr 
Josephson junction and we compare the results. Finally, we will discuss further open questions 
in a Conclusion. 

2. Fluxons and Semifluxons in Superconductivity 

The Josephson effect is a well established phenomenon in the solid state physics. A Josephson 
junction (JJ) consists of two weakly coupled superconducting condensates. JJs are usually 
fabricated artificially using low- or high-Tc superconducting electrodes separated by a thin 
insulating (tunnel), normal metal, or some other (exotic) barrier. JJs can also be present 
intrinsically in a anisotropic layered high-Tc superconductors such as Bi9Sr9Cai Cu90« [l32l 

123. 

The dc Josephson effect (flow of current through a Josephson junction without producing 
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voltage drop, i.e. without dissipation) is expressed using the first Josephson relation, which in 
the simplest case has the form 

= Jc sin (0), (1) 

where Ig is the supercurrent flowing through the junction, is the critical current, i.e. the 
maximum supercurrent which can pass through the JJ, and (p = 62 — Oi h the difference 
between the phases of the quantum mechanical macroscopic wave functions 2 = y^e'^^-^ 
of the superconducting condensates in the electrodes. 

Recent advances in physics and technology allow to fabricate and study the so-called 
TT-Josephson junctions — junctions which formally have negative critical current < 0. 
This can be achieved by using a ferromagnetic barrier, i.e. in Superconductor-Ferromagnet- 
Superconductor (SFS) jMl ESI [36l [37l |38l or Superconductor-Insulator-Ferromagnet- 
Superconductor (SIFS) [[39l l40l structures. One can also achieve the same effect using a 
barrier which effectively flips the spin of a tunneling electron, e.g. when the barrier is made 
of a ferromagnetic insulator ||4TI . of a carbon nanotube ll42l or of a quantum dot created by 
gating a semiconducting nanowire ll43ll . 

The change in the sign of a critical current has far going consequences. For example, 
analyze the Josephson energy (potential energy related to the supercurrent flow). In a 
conventional JJ with /c > 

U{ct>)=Ej{l-cos(t>) (2) 

and has a minimum at = + 2?™ (the ground state), where Ej = $o/c/27r is the Josephson 
energy. If Ic < 0, we define Ej = <l>o|/c|/27r > and 

U{(j))=Ej{l + cos<P). (3) 

Obviously, the minimum of energy is reached for = tt + 27rn. Thus, in the ground state (the 
JJ is not connected to a current source, no current flows through it) the phase drop across a 
conventional JJ with Jc > is = + 27m, while for a junction with < it is = tt + 27m. 
Therefore, one speaks about "0-JJs" and "tt-JJs". 

Further, connecting the two superconducting electrodes of a vr-JJ by a not very small 
inductor L (superconducting wire), the supercurrent oc n/L will start circulating in the loop. 
Note, that this supercurrent is spontaneous, i.e., it appears by itself, and has a direction 
randomly chosen between clockwise and counterclockwise [44] . The magnetic flux created 
by this supercurrent inside the loop is equal to $o/2, where $0 = h/2e ~ 2.07 ■ 10^^^ Wb 
is the magnetic flux quantum. Thus, the vr-JJ works as a phase battery. This phase battery 
will work as described, supplying a supercurrent through the loop with inductor, provided the 
inductance L ^ ^o/h- If the inductance L is not that large, the battery will be over loaded, 
providing a smaller phase drop and supporting smaller current. For very small inductance the 
battery will stop working completely. 

Similar effects can be observed in a tt dc superconducting quantum interference device 
(SQUID: one 0-JJ, one vr-JJ and an inductor L connected in series and closed in a loop) or in 
0-vr Josephson junction. Let us focus on a latter case. 
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Figure 1. Sketch of a O-tt JJ and profiles of (a) the phase (t>{x), (b) the magnetic field 
oc d(t){x)/dx and (c) the supercurrent density js{x) oc sin0(x) corresponding to a semifluxon 
(black) and antisemifluxon (gray). 



Consider a long (along x) Josephson junction (LJJ) one half of which at x < has 
the properties of a 0-JJ (critical current density jc > 0) and another half at x > has the 
properties of a vr-JJ (critical current density jc < 0). Long means that the length is much 
larger than the so called Josephson length Aj, which characterizes the size of a Josephson 
vortex; typically Aj ~ 10 — 20 fim. What will be the ground state of such a O-tt LJJ? 
It turns out that if the junction is long enough (formally infinitely long), then far away 
from the O-tt boundary situated at x = 0, i.e., at x ^ ±oo the phase (p will have the 
values or ivr (we omit 27m here), while in the vicinity of O-tt boundary the phase (p{x) 
smoothly changes from 0(— cx)) = to 0(+oo) = ±7r, see Fig. d^. The exact profile can 
be derived analytically [|45ll46l[3n . Since the phase bends, the local magnetic magnetic field 
H oc d(f)/dx will be localized in the vicinity of the O-tt boundary and carry the total flux 
equal to ±$o/2, see Fig.[ll). The sign depends on whether the phase bends from (^(— oo) = 
to 0(+oo) = +7r or to (j){+oo) = —n. Thus, such an object is called a semifluxon or an 
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antisemifluxon. If one analyzes the Josephson supercurrent density flowing though the barrier 
js{x) = jc{x) sin</)(x), one can see in Fig. [l]: that the supercurrent has different directions 
on different sides from the O-tt boundary. Since we do not apply any external current, the 
flow of current should close in the top and bottom electrodes, i.e. the supercurrent circulates 
(counter)clockwise in the case of (anti)semifluxon. Thus, a semifluxon is a Josephson vortex 
of supercurrent. It is pinned at the O-tt boundary and has two degenerate ground states with 
the localized magnetic field carrying the flux ±$o/2. 

Semifluxons in various types of JJ has been actively investigated during the last years. 
In fact, the first experiments became possible because of deeper understanding the symmetry 
of the superconducting order parameter in cuprate superconductors. This order parameter 
with the so-called ci-wave symmetry is realized in anisotropic superconductors, such as 
YBa2Cu307 or Nd2_a;Ce^.Cu04. It allowed to fabricate O-vr grain boundary LJJs [El SHI 
and, later, more controllable d-wave/s-wave ramp zigzag JJs ll49l l50l and directly see and 
manipulate semifluxons using a SQUID microscope [[5T1I521 . 

Semifluxons are very interesting non-linear objects: they can form a variety of ground 
states [I53l[5l[55l[56l, may flip emitting a fluxon [[5711551 EH, or be rearranged W\ 

by a bias current. Huge arrays of semifluxons were realized [[STI and predicted to behave as 
tunable photonic crystals [[60l . Semifluxons are also promising candidates for storage devices 
in the classical or quantum domain and can be used to build qubits [[611 as they behave like 
macroscopic spin 1/2 particles. 

Now, an interesting question arises: Can one realize vr or even O-tt JJs in an atomic 
BEC? In the latter case, the degenerate ground state corresponding to a semifluxon should 
have a non-trivial spatial phase profile and semifluxon physics can also be studied using BEC 
implementation. 

3. Semifluxons in Bose-Einstein Condensates 

Here, we will address this question and examine a configuration were the two-state atoms 
are trapped in a quasi-one dimensional cigar-shaped trap with an additional superimposed 
double well potential in the longitudinal direction. The spatial localization of the two-state 
atoms inside the double well potential leads to two internal atomic Josephson junctions that 
are driven via an optical, complex valued "O-tt" laser field and they are motionally connected 
via tunneling. 

First, we will describe details of the model and introduce the Hamiltonian of the system. 
Then, we will examine the classical limit of the field theory and study the ground state of 
the Gross-Pitaevskii equation. Finally, we will exploit the fact that the spatial wavefunctions 
are localized inside a deep double well and study a simple four-mode quantum model derived 
from a Wannier basis state representation. 
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3.1. O-TT-junction in a Bose-Einstein condensate 

To model the O-vr-Junction ll29l [30l in a BEC, we are guided by the the condensed 
matter physics setup depicted in Fig. [2l As in the single atom case, we replace the two 
superconductors by an atomic two-level BEC in a cigar shaped trap. The two states of the 
atom, i.e., the excited state |e) and ground state \g) couple via a position dependent Rabi 
frequency Vto{x), which exhibits a phase jump at the origin of the x-axes 




a) 






Figure 2. Analogy of a O-tt Josephson junction in a superconductor (a) with an atomic two- 
level BEC (b) in a double well trap V{x) with a position dependent Rabi frequency ilo{x). 
The excited atomic level |e) is separated from the ground level \g) by the detuning 26. 

In this quasi one-dimensional scenario, we will represent the two-state atoms by a 
spinorial bosonic quantum field ^, which satisfies the commutator relation 

= M(x-y). (5) 

This field can be decomposed in any complete single particle basis \a,l), which resolves the 
spatial extent of the field and the internal structure of the atoms, i.e., 



a={e,g} 1=0 

and we denote the corresponding discrete bosonic field amplitudes by ao-,i- Here, a 
characterizes the internal states by (e, g) and the external motion in the double well potential 
by a quantum label /. The dynamical evolution of the atomic field is governed by following 
Hamiltonian 



H 



< + V{x) + 



5 

ni{x) 



fio(x) 
-5 



^!(x) 



(7) 



.X 



In here, we use dimensionless units, in particular we have set /i = 1 and the mass of the atom 
m = 1. The energy consists of the single particle energy in a trap V{x), which is identical 
for both species, the electric dipole interaction of the two- state atom ll62l . as well as a generic 
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collision energy proportional to the coupling constant g = gee = ggg- To simplify the analysis, 
we have deliberately set the cross-component scattering length geg = 0- No unaccounted loss 
channels are present. Therefore, we have number conservation 

H,n]=0, (8) 

as a symmetry. If we denote a generic state of the many -particle system in Fock-space by 

oo 

I^W) = 5^V^n(t)|n=(no,ni,...)), (9) 

n=0 

then one can obtain the dynamics of the system most generally from the Lagrangian 
formulation 

c[Mt),Mt)] = mrndt-Him)- m 

In this field theory, the canonical momentum is given by vTn = = iifj'l^. From the Hamilton 
equation tx^ = —j^, one recovers the conventional Schrodinger equation in Fock space 

idt\^) = H\^). (11) 

The Lagrangian approach is obviously a central concept in the path integral formulation of 
quantum mechanics [l63ll . However, it is also of great utility in the approximate description 
of the dynamics if we connect it with concepts of classical mechanics as we will see in the 
following. 



3.2. Spatially extended classical model: the Gross-Pitaevskii equation 

The classical limit of the field equations [|2l can be recovered quickly by approximating the 
state of the system by a coherent state 

i!,{x)\iP)=^,{x)\iP). (12) 



Within this approximation, we obtain from the Lagrangian of Eq. (ITOI) the two component 
Gross-Pitaevskii equation = {'ipe{x,t),'ipg{x,t))~^ 



-dl + V{x) 



VLq{x) 

mix) -5 + 2g\i,g\'' 



5 + 2g\ilj, 

^0 



(13) 



For a macroscopically occupied field this equation models the spatial evolution of the coupled 
Josephson junctions very well |[l7l[Il[BliQl|2lll22ll23ll2ll3ni. 



3.3. Discrete quantum model: two coupled Josephson junctions 

To gain more insight into the quantum properties of the groundstate of the system ||6TI . one 
can decompose the field into its principal components and disregard small corrections. A 
double well potential can be considered as the limiting case of a periodic lattice. While even 
and odd parity modes relate to delocalized Bloch-states in a periodic system, left and right 
localized modes i.e., (po and right (pi, resemble the Wannier basis. These localized modes are 
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depicted in Fig.[3l With respect to these basis states, we can approximate the field with four 
modes 



(14) 
(15) 



^000 (a;) +^i0i(x) +5i!g. 

The four bosonic amplitudes {^q, Cq, 5?]^, e^^} satisfy the usual commutation relations and we 
will disregard the small corrections of order 5'^„ 




Figure 3. From even and odd parity modes of a double-well potential, one can construct the 
left 00 ™d right 0i localized Wannier states of the system. 



If this approximate field is substituted in the Hamiltonian, Eq. (|7]), we can exploit the 
orthogonality of the wave functions and obtain the two-body matrix elements (pij^i 

/•oo /»oo 

6ij= / dx0i(x)0j(x), (t)ijki = / dx0i(x)0j(x)0fc(x)0/(x). (16) 



Out of the sixteen combinations for (pijku we only retain the physically most relevant 
contributions and disregard others deliberately. This leads to the following model Hamiltonian 
for two coupled Josephson junctions 

H = k{ele^ + e\e^+glg^ + g\g^) (17) 

+ ^(4^0 - ^0^0) - qIqi) 

+ fio(4^o + ^Jeo) + ^Mdi + 9\ei) (18) 
+ ai^lele^e^ + glgl%% + e\e\e^e^ + glglg^gi) 

where all coupling constants are implicitly rescaled by the corresponding single particle or 
two-body matrix elements and the new parameter A measures the spatial hopping rate between 
the sites. 
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Figure 4. Coupling pattern of two coupled Josephson junctions JJq and JJi with, fii = 
— r^o- The classical phase Q refers to a particle current between the spatial sites, while the 
phases 0o and 0i are currents within the two-state atoms. 



SA. Fock-space representation of the four mode model 

In principle, it is possible to solve the four-mode Schrodinger equation in Fock space by 
projecting it on the A^-particle sector 

oo 

1^, iV) = 5(iV - iV) ^ ^„|n = {nl, nl n\, n\)). (19) 

n=0 

The number constraint on the state remains valid throughout the time-evolution as number 
conservation is encoded into our Hamiltonian from the beginning in Eq. ([8]). This reduces 
the discrete d = 4 dimensional eigenvalue problem to an effective three-dimensional problem 
with nontrivial boundaries. We have illustrated the finite support of the amplitude field in 
Fig.lH It is a (d — 1) -dimensional simplex embedded into a rf-dimensional Fock space. The 
full analysis of this problem is an interesting problem in its own right and will be presented in 
a forthcoming publication. 

3.5. The classical limit of the four mode model 

It is not necessary to solve the four-mode problem in Fock-space to understand the principal 
features of the equilibrium configuration. Thus, we will again resort to the classical 
approximation and use the number-symmetry broken coherent state approximation for the 
quantum state 

1^) = |q; = (eo,5fo,ei,5fi)'^), ei\ii) = ei\ii) , gi\'ilj) = gi\ilj) . (20) 
The dynamics is simply obtained from the Lagrangian 

i d 

C{oL, a) = {a\idt — H\ci) = ia*a — H{a, a*) — 2'^'^^'^' ^*^' ^^^^ 

if we introduce the classical Hamilton functions 'H{a.,a*) = {a.\H\a.) and number 
expectation values A/'(q;, a*) = {oc\N\a), as 

n{a, a*) = A (e*ei + ejeo + g^gi + glgo) (22) 
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Figure 5. A three dimensional simplex represents the four-mode Fock space for = 32 
particles. The axes are labeled in a generic lexicographical order (no, ni, ^2) and implicitly 

= N — tiQ — ni — 71,2. 



+ 5o(ieop-kr)+5i(ieir-kn 



^(|eo 

|2 



^ol 

2 



+ |ei| 

|2 



eol + Ifi'ol + |ei 



+ 1^71 1 



If the dynamical coordinate is a, then we find the canonical momentum as 



TTfc 



dak 



la 



Consequently, variables and momenta obey the Poisson bracket {aj, tt^} 
construction, number conservation is satisfied dynamically as 



dt 



{HM] = 



and the Hamilton equations of motion for the coordinates read 

d 



dt 



-a 



( 5o + 2^|aoP 

A 




ilCoL 



'jk- 



^0 

f 2(7|ai|2 

A 



Si 





A 

-5i + 2^|a3|2 J 



(23) 

(24) 
By 

(25) 
(26) 



A 

-<5n + 2(7laiP 

'2^|a2p 

In classical mechanics, we can deliberately choose new coordinates and momenta to 
accounts for symmetries of the Hamiltonian. If a new variable of a canonical transformation 
matches a conserved quantity, it follows that the conjugate variable becomes cyclic. In this 
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spirit, we will introduce the following pairs of action-angle variables: ($,Ar) measures 
the global phase and the total particle number M = Uq + Uq + + of the system, 
{9, Ai) measures the relative phase between left and right sites and the population imbalance 
Ai = nQ+nff — {nl+nl) inbetween, (0o,-^o = ''^o~''^o)' ^^'^ = nl—nf) measure the 

relative internal phase of the atoms on each site and the corresponding population difference. 
This coupling scheme for the phases and population imbalances has been illustrated in Fig.lH 
By inverting the population relations, one finds the individual occupation number nf per site 
as 

nl = \{^-M + 2Mi), nf = l{^-M-2Mi), ^ ' 

Finally, we can use these physical coordinates in a canonical transformation from complex 
amplitudes to real action-angle variables 

g^^g-i(f+e+<^o) /^^ g^^^-i{^+e-M^^ (28) 



For later use it is also useful to introduce the auxiliary phases 6*^, which are global phases of 
the subsystem on site i: 

9^ = ^ + 6, ei = ^-e. (29) 

By substituting field amplitudes into the Lagrangian of Eq. (|2T1) . one can again identify the 
variables and corresponding canonical momenta as 

dC dC . , dC . , dC . , 

^=M, ^ = M, ^ = Mo, ^ = Mi. (30) 

With this new coordinates and momenta we obtain a Hamilton function 

n' = n'{e, 0o, <j)i,^,M,Mo,Mi) Od 

= ^ cos {26 + 00 - 0i) A/(Ar + M + 2Mo) {AT - M + 2Mi) 
+ ^cos{2e-(f)o + (j)l)^/{^f + M-2Mo){^^-M-2Ml) 



fio cos (20o) J (AT + 7W) 2/4 -MI + 6oMo + S^Mi 



+ fii cos (20i) ^J{^^-M)yA-Ml + ^{^^^ + M^ + 2mI + 2mT) . 

For the parameters 6 = 1, = ^Mi = -1, A = 0.1, g = 0.005, A/" = 100, we have 
numerically minimized the energy of Eq. (|3T1) and find two minima E = Eq, Ei. Please note 
that the amplitudes of the solutions a. are real valued. Due to the invariance of the solutions 
under global phase change, we can deliberately modify them and choose the global phase of 
the left site as the common ground 6^0 = 0. In Tab. [H we have listed the energies, the real 
valued amplitudes and the phases calculated according Eqs. (l28l) and (l29l) . 
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Table 1. Two stationary solutions of the Hamilton equation, energies E^, Ei, field amplitudes 
a. = (eo, (?0j ei, gi) and phases. The global phase of the left site, i.e., 6q = was chosen as 
the common ground. 
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Figure 6. Comparing a spatially extended 0-7r-junction in a superconductor (left panel) 
calculated according to the sine-Gordon equation IMl to an analog, two site junction in an 
atomic BEC (right panel). We have plotted the scaled global phase 9i vs. the site index 
i = 0,1. 



4. Conclusion and Perspectives 

In the present article we have briefly summarized the status of the fluxon and semifluxon 
physics in superconductivity. In particular, we were focusing on an effect that exists in a 
long O-TT geometry where two Josephson junctions were in the groundstate and one obtains a 
Josephson vortex of fractional magnetic flux, pinned at the O-tt boundary. Depending on the 
preparation procedure, this "classical" state of the superconducting device can occur in two 
different configurations: with the magnetic flux equal to +$o/2 or — $o/2. 

This arrangement of superconducting junctions has been analyzed and transfered to 
the context of an atomic bosonic quantum gas, where two-state atoms in a double well 
trap are coupled in a similar fashion. There, the optical O-vr junction is represented by 
the left/right localized internal atomic Josephson junctions and a jumping phase of the 
complex-valued Rabi-frequency. We have derived a simple four-mode model for this case 
and showed that in the "classical" approximation it qualitatively resembles the semifluxons 
seen in superconductivity. 

In the superconducting case, one observes a smooth spatial behaviour of the phase across 
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the O-vr boundary. In the atomic case, this will emerge also by a more realistic modeling for the 
Rabi-frequency and spatial motion. This is currently under investigation. Eventually, the four- 
mode, O-TT Josephson model will be also instrumental in examining the quantum properties of 
these macroscopic semifluxon states and we will explore the macroscopic tunneling between 
them lIMI . This is also work in progress and will be reported in a forthcoming publication. 
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